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Two identical chaotic oscillators that are mutually coupled via time delayed signals show very 
complex patterns of completely synchronized dynamics including stationary states and periodic as 
well as chaotic oscillations. We have experimentally observed these synchronized states in delay- 
coupled electronic circuits and have analyzed their stability by numerical simulations and analytical 
calculations. We found that the conditions for longitudinal and transversal stability largely exclude 
each other and prevent e.g. the synchronization of Pyragas-controlled orbits. Most striking is the 
observation of complete chaotic synchronization for large delay times, which should not be allowed 
in the given coupling scheme on the background of the actual paradigm. 
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I. INTRODUCTION 

Chaos control and chaotic synchronization have been 
one of the most attractive fields of applied nonlinear 
dynamics for almost two decades in physics P, 
engineering Ql, chemical and biological systems 
@ and others. Various strategies aiming at the 
suppression of undesired chaotic oscillation or at the 
control of unstable periodic orbits have been developed 
and have meanwhile been well understood. For a review 
see Q- The importance of chaotic synchronization in 
complex systems has only gradually been clarified [8l-[TT| . 

Chaotic synchronization means that two or more 
dynamical systems, which are coupled with each other, 
can correlate their dynamic behavior, even in the case 
when each system shows irregular behavior. In general, 
various grades of synchronization can be obtained 
depending on the coupled elements but also essentially 
on the form and strength of the couplings. Here we will 
focus on the strongest form of synchronization called 
complete or identical synchronization, which means that 
the trajectories of the coupled systems coincide exactly. 
In the case of two identical chaotic systems complete 
synchronization can be achieved by non-invasive cou- 
plings. This means that all coupling forces vanish as 
soon as the synchronized state has been reached, and 
the synchronized state will be in exact agreement with 
the previously uncoupled states. 

Recently the importance of latency effects in the 
coupling mechanism have attracted major attention 
both from the aspect of practical applications and for 
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principle reasons. For example, when developing a high 
speed cryptographic communication system based on 
chaotic synchronization [l^ . it will be most important 
to understand the influence of latency times within 
the whole complex system. Many recent investigations 
of this problem refer to coupled lasers, e.g. |13l - {2g |. 
because of their fast characteristic time scales, but also 
owing to the intrinsic applications of delay lines to make 
a single laser element chaotic. 

Delay effects may, however, substantially change the 
mechanism of chaotic synchronization. It is evident that 
for two irregular time-dependent signals x(i) and y{t) 
any time delay in the coupling forces has the tendency 
to create a mismatch in time and will impede the 
state of complete synchronization. This was the reason 
for the previous opinion that chaotic systems with 
mutual time-delayed couplings could never synchronize 
completely. Obviously they cannot synchronize in a 
non-invasive way, like unstable periodic orbits [2]| or like 
chaotic systems with instantaneous @ or unidirectional 
[2^ couplings. Instead non-zero delay in bidirectional 
couplings will always result in strong control forces which 
have to overcome the deviations of identical irregular 
signals with some time lag. This makes the coupling 
invasive. In this context it has been recently realized 
that more complicated coupling schemes including 
additional self-feedback may partly remove the influence 
of delay from the eqs. of motion of the coupled systems 
[23| , but even then the principle conflict between delayed 
couplings and the synchronization of irregular signals 
is still maintained. Of course there will be further 
constraints, and we may expect that the tendency of 
nearby trajectories to exponentially diverge cannot be 
suppressed by delayed couplings for very large delay 
times [H]. 
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FIG. 1. Coupling scheme. 




In this paper we focus on the most simple coupling 
scheme and investigate two identical low-dimensional 
chaotic oscillators that are mutually coupled via a single 
time-delayed component. We study, both experimentally 
and analytically, the mutual and auto-synchronization ef- 
fects that appear in this configuration. Our main inter- 
est will be directed to the following problems: (i) Can 
such chaotic oscillators completely synchronize at non- 
zero time delay, and what are the limits for delay time? 
(ii) Is it possible to control unstable periodic orbits em- 
bedded in these chaotic oscillators and to synchronize 
them simultaneously? 



II. EXPERIMENT 

For our experimental investigations we designed two 
identical autonomous electronic oscillators from analog 
components including active elements. A diagrammatic 
view of our setup is shown in fig. [21 

Following Kirchhoff 's rules the equation of motion of 
this circuit can be written as 



C/i ^ ujo {-a U1-U2- C/3) 
U2 = Wo {Ui + {a-a) U2) 
U3=I3Id{Ui,U3)^iU3 . 



(1) 



The typical time scale of chaotic oscillation is determined 
by the parameter luq which depends on the RC values 
of integrating op-amps and - for practical reasons - was 
chosen to be in the order of 1ms. The dimensionless 
parameters a and a are determined by the ratios of 
resistors. For details see [2^. Here, a is considered as 
the control parameter of the oscillator. The nonlinear 
characteristics of the diode D (DC 1N4007) provides the 
nonlinearity of the circuit. The dynamics of each single 
oscillator reflects almost perfectly that of a Rossler 
oscillator, in spite of having replaced the multiplicative 
nonlinearity of the original Rossler model by the nonlin- 
ear characteristics of a diode. Circuits with such a type 
of nonlinearity have turned out [25| to be more robust 



and less sensitive to drifting offsets than circuits includ- 
ing electronic multipliers. Finally, a slight additional 
damping term a = 0.05, not present in the original 
Rossler model, has been attributed to the Ui and U2 
components in order to avoid a drift in the integrator 
output signal due to some slight unavoidable offset at the 
input. The corresponding electronic components of both 
'identical' oscillators were compared and carefully se- 
lected to guarantee that any deviations remain below 1%. 




FIG. 2. Design of a single Rossler-type electronic oscilla- 
tor. Values of components: Ri = 2.7kfl, R2 = 3.6k^l, Ri = 
7.5fcn, R4 = lOkQ, R5 = 13kn, Re = 15kQ, R7 = SSfcfi, Rs = 
200kQ. Ra is variable — 50^57 and determines the system 
parameter a according to Rg — lOkil/a (accuracy of resis- 
tors 1%). Capacitators: C — lOnF (accuracy 5%). Type of 
diode: 1N4007 from DC Components. Type of operational 
amplifiers: TL084. 



The mutual delayed couplings between the oscillators 
were designed in the simplest possible way via a single 
component. We chose the U2 components of both oscil- 
lators (labeled x and y). The input signals C/2,2; and U2,y 
were sent through separate digital delay lines and multi- 
plied by a gain factor k to compose the coupling signals 
k {U2,y[t - t) - U2At)) and k (?72,.(t - r) - U2,y{t)), 
respectively. In order to obtain well-defined delay sig- 
nals the feedback loops generating the control signal are 
switched on at about one cycle later than the oscillator. 
Measuring time in units of lo^^ and all voltages in volts 
the coupled system is described by the normalized equa- 
tion of motion 

Xi = —axi — X2 — X3 

±2 ^ xi + {a - a) X2 + k (?/2,t - X2) 

X3 = gixi,X3) - CX3 



(2) 



yi = -ayi - 2/2 - 2/3 

2)2 = yi + (a - a) y2 + k {x2,t - 2/2) 

2)3 = 5(2/1,2/3) - C2/3 . 
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FIG. 3. Diode characteristics in the network of resistors 
Ri, Rq, Rj. In order to simulate the situation in the running 
circuit, a sinusoidal test voltage UM{t) with 5V amplitude 
and a frequency of 1.6kHz was applied corresponding to the 
dynamic input from Ui and U3. At the transition point be- 
tween locking and conduction {Um ~ 3V) a small undershoot 
occurs due to the diode capacity, which is negligible for the 
dynamics of the circuit. 



Here, Xi and yt denote the normalized components 
Ui.x/iV and Ui,y/lV, respectively. The control param- 
eter was fixed to a = 0.265, a = 0.05, and c = 7/cjo = 
2.82. The normalized diode characteristic g{-), which 
induces the folding process of the Rossler-type dynam- 
ics, was approximated by a piecewise linear function 
5(2:1, Z3) = f3 {\z\ + z) with z ~ zi + Z3/2 - zthr- f3 and 
Zthr were obtained from a fit of this function to the exper- 
imentally obtained profile shown in fig. [3] The setting of 
system parameters, delay time and coupling amplitudes, 
and the acquisition of data are handled by a separate 
module, allowing to run the whole experiment automat- 
ically. 



III. OBSERVED PHENOMENA 

In general there is a large variety of dynamical phenom- 
ena arising from this simple setup due to the interplay of 
nonlinearity and delay. In order to get an overview, we 
first study the dynamical states occurring on variation of 
coupling strength and delay time. Figure 0] shows an ex- 
perimental diagram of different dynamical states. Only 
those parameter combinations are marked for which com- 
pletely synchronized dynamics is observed, and different 
colors correspond to different dynamic states, e.g. syn- 
chronized orbits of various periodicities as well as syn- 
chronized chaos, and - as a trivial limit - also 'synchro- 
nized' (i.e. identical) fixed points of both systems. The 
diagram shows that for the parameter range investigated 
complete synchronization is rather an exception than the 
rule. 



3 




0.1 0.2 0.3 0.4 

k 



FIG. 4. (Color online) Experimental diagram of synchronized 
states. NS: Not completely synchronized, IM: intermittent, 
AP: aperiodic, FP: fixed point, P: complete periodic synchro- 
nization. 



For T ~ 0, which represents the conventional case with- 
out delay, we have complete chaotic synchronization 
when the coupling strength k exceeds a critical value kc- 
For delay times slightly larger than zero, we still observe 
complete synchronization but apart from synchronized 
chaotic states we also observe synchronized orbits and 
fixed points. It was counterintuitive to see that peri- 
odic orbits do not or almost not synchronize for r = T 
or integer multiples of T, but slightly aside there are 
stripes of complete periodic synchronization. There is 
also no complete synchronization for half-integer mul- 
tiples of T. The most remarkable phenomenon of CS 
is a boomerang-shaped domain of complete chaotic syn- 
chronization around r — 3T. The following sections are 
devoted to analyze these patterns of complete synchro- 
nization. 



IV. THEORETICAL APPROACH 

In order to discuss the problem in more general terms 
we start from two identical oscillators x = f(x) and 
y = f(y), which are symmetrically coupled by a bidi- 
rectional interaction which is proportional to the instan- 
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tancous difference of their trajectories: 

x-f(x) + K(y-x) 
y = f (y) + K(x - y) . 



(3) 



This conventional type of coupling can be considered as 
the limiting case of zero time delay. The coupling acts as 
an additional dissipation, and it is well-known that there 
exist coupling matrices K such that the completely syn- 
chronous state x(i) = y{t) becomes stable. If the signal 
arriving from the other oscillator is delayed in time, e.g. 
due to finite propagation velocity in a long transmission 
line, then the ordinary differential cqs. ^ turn into a 
system of time-delayed differential cqs. 



X = f (x) + K(y, - x) 
y = f (y) + K(x, - y) 



(4) 



where Xt- = x(t — t) and yr = y(i — t). Although such 
a system is very difficult to tackle, we sec immediately 
that the completely synchronous state x{t) = y{t) still 
represents a solution of cqs. Q. Intuitively it is not clear 
whether such a state can be stabilized in presence of the 
delay terms. As we already discussed above, there is a 
principle difference between the synchronous solution of 
eq. ([3]) and that of eq. (|4]). In the case of instantaneous 
interactions the coupling forces vanish as soon as the 
completely synchronized state has been reached, i.e., 
this kind of coupling is non-invasive. Time-delayed cou- 
plings, however, are generally invasive. When complete 
synchronization has been obtained, the time-delayed 
signal of the remote system cannot be distinguished from 
the local signal delayed in time. Thus we may consider 
the coupling as a virtual self-feedback, and the resulting 
dynamics of the coupled systems eqs . Q represents that 
of a single delay-system (see e.g. |26|). 

Stability analysis of synchronized solutions generally 
starts from the transformation into appropriate normal 
coordinates. According to the symmetry of our problem 
we define the longitudinal coordinates u = (x + y)/2 
and the transversal coordinates v = (x — y)/2. In 
terms of theses coordinates the synchronization mani- 
fold is determined by u = x = y and v = 0. For 
small synchronization errors we use the linearized form 
f (u ± v) = f (u) ± Df (u)v, so that up to first order in v 
the eqs. ^ can be written as 



u = f(u) +K(u^ - u) 

V = Df (u)v + K(-v^ - v) 



(5) 



The equation of motion for u just reflects that of a single 
original oscillator with time-delayed self-feedback, which 
is in exact analogy with the Pyragas control scheme. 
This form again underlines the - in principle - invasive 
character of the present coupling due to the additional 
delay term occurring in the equation. For zero delay, 
however, the additional control-induced term vanishes 
and the coupling type becomes non-invasive. The 



transversal coordinate, v, builds up a special tangent 
linear system (TLS) driven by u. The stability of its 
solution V = is determined by the interplay between 
chaotic drive and coupling-induced damping terms, 
which include both instantaneous and time-delayed 
values of v. 

The stability properties of both eqs. depend in a very 
complex way on both the delay time r and the strength 
and type of the coupling described by matrix K. The 
basic criterion for complete synchronization is the stabil- 
ity of the transversal solution v = 0. If the solution of 
the longitudinal eq. has at least one positive Lyapunov 
exponent we will observe the synchronization of chaotic 
trajectories. However, the presence of the self-feedback 
term may also result in stable periodic solutions or even 
fixed points, when the largest Lyapunov exponent be- 
comes zero or negative, respectively Yet, since both the 
longitudinal and the transversal eq. depend on the same 
coupling parameters K and r we have to explore whether 
the conditions for such regular solutions for u{t) are con- 
sistent with the stability of v = 0. In the following sec- 
tions we will analyze these different types of solutions 
separately. 



V. STABILITY OF THE FIXED POINT 

A free running Rossler-type system like that of cq. ([T]) 
has an unstable fixed point close to the origin u = 0. 
Its exact position depends on the nonlinearity .g(x). In 
our experimental system the fixed point coincides exactly 
with the origin: g{0) = 0. We are interested in the 
longitudinal and transversal stability properties of the 
synchronized fixed point, i.e. of the solution where the 
two oscillators lock in their steady states due to the cou- 
pling. Note that in the case of steady states our coupling 
scheme will be non-invasive, since the delay term in the 
longitudinal part of eqs. ([S]) vanishes. Introducing very 
small deviations du and Sv from this common fixed point 
in both the longitudinal and transversal components its 
stability properties can be obtained through the following 
linearized equations: 



Su = Df (0)(5u -I- K{+Sur - Su) 
6v = Df (0)(5v + K{-Svr - Sv) , 



(6) 



where Df(0) denotes the Jacobian matrix evaluated at 
the origin. It creates the common drive for both longitu- 
dinal and transversal direction. Due to the simplicity of 
stationary states the linearized equations decouple, and 
we have two independent tangent linear systems that de- 
termine the stability of the synchronous fixed point. The 
Jacobian of our system reads explicitly 



Df(0) 




(7) 
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FIG. 5. Longitudinal (/ / /) and transversal (\\\) stability 
domains of the fixed point. 



To further simplify our analysis we neglect the Su^ and 
Sv3 components, because after short transients both will 
rapidly relax to zero: Su^ oc 6V3 cx exp(— ct). Only the 
dominant oscillating part remains relevant for stability 
and will be treated in terms of a complex- valued normal 
form: 



i = XqZ + K{ + Zr — z) 
W ~ XqW + K{ — Wr — w) 



(8) 



where z = 5ui + iSu2 denotes the longitudinal and w = 
6vi + iSv2 the transversal linearized coordinate. Aq = 
(a/2 ~ a) + i is one of the complex conjugate eigenvalues 
of the remaining 2x2 Jacobian matrix, and the coupling 
is At = fc/2. These results are approximations of the real 
dynamics for the case of strong mixing of the components 
due to rotation. An exact transformation is not possible 
because Df and K do not commute. Our approximation 
is valid as long as the coupling strength is not too large. 
The corresponding eigenvalues of eqs. (|5]) are determined 

by 



A|| =Ao-K(l-e-^ii") 



A^ = Aq — k{1 + e 



-A, 



(9) 



The synchronized fixed point is stable if all 5R(A) < 0. 
Figure [5] shows the stability domains in the fc, r-plane. 
The domains of transversal stability arc centered at in- 
teger multiples of T, while the domains of longitudinal 
stability arc centered at half-integer multiples. There is 
only a very small overlap between both kinds of domains 
at a quarter-period delay t = T/4. In fig. [5] we have il- 
lustrated the fc-dependence of both the longitudinal and 
transversal Lyapunov exponents at a half-period delay 
T = T/2. For the longitudinal component we recognize 
the profile familiar from the Pyragas control of a peri- 
odic orbit with 7r-torsion [27[. But for the transversal 



FIG. 6. Real part of the largest eigenvalue 5R(A) from eq. (|8]) 
for both longitudinal and transversal subspace at a half-period 
delay r = T/2. 



component the delay term is of opposite sign, so it coun- 
teracts with the instantaneous term - just like in the case 
of torsion-free orbits, which cannot be stabilized by ordi- 
nary Pyragas control. 

The overlap at the quarter-period delay arises from 
a 7r/2 phase rotation, for which the opposite signs in 
eqs. (O lead to complex conjugate eigenvalues. We 
found that such an overlap can only occur as long as 
3^(Ao)t < 0.326 [2^. This explains the missing over- 
laps at larger delay times. Altogether, the experimental 
observation of synchronous fixed points in a very small 
regime at r « T/4 only is well understood in terms of 
our stability analysis. 



VI. PERIODIC SYNCHRONIZATION 

Each chaotic oscillator has included an infinite set of 
unstable periodic orbits. In this section we focus on the 
low-periodic orbits, which are embedded in the chaotic 
attractor of an uncoupled oscillator. There exists at least 
one period-1, period-2, etc. orbit with a period, which 
is close to an integer multiple of the mean cycle time 
T « 27r/wo. Moreover, since in our circuits these unsta- 
ble orbits arise from a period-doubling route to chaos, 
they all have a finite torsion, which means that their 
neighborhood undergoes a 7r-flip during each cycle r„. 
Therefore these orbits should be accessible to Pyragas 
control under the common restrictions. In fact, when we 
applied Pyragas control with r = T„ and appropriate 
coupling strength to a single oscillator, the stabilization 
of such orbits could be observed in experiment. But when 
trying to synchronize two identical Pyragas-controllcd 
orbits by means of our mutual delay-coupling scheme 
eq. ([2) - which should be non-invasive in this case, too 
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- wc always failed. This becomes evident from fig. |4] for 
T = T, 2T, and AT. Note that the narrow synchroniza- 
tion windows observed at smaller fc-values merely refer 
to delay-induced periodic solutions for which the feed- 
back force does not vanish. So we conclude that for our 
coupling scheme the non-invasive synchronization of sta- 
bilized periodic orbits is not possible. This can be ex- 
plained by analogy with the synchronization mechanism 
of fixed points. Let denote the unstable periodic or- 
bit of an uncoupled oscillator with ^{t + T) = ^{t). Then 
for T = nT, n integer, the trajectory x(i) = y{t) = ^{t) 
is a solution of eqs. ^ which we want to analyze. In- 
troducing small longitudinal and transversal deviations 
5u and 6v from the orbit, we arrive at almost the same 
two tangent linear systems as in eqs. ^ but with a time- 
periodic drive instead of a constant one. 

Su = Bi(^(t))6u + K(+dUr - 6u) 

(10) 

6v = Df (4(t))Jv + K{-5vr - 6v) . 

The longitudinal equation represents exactly the Floquet 
problem well-known for Pyragas control. So within the 
synchronization manifold the trajectory shows the same 
stability properties as an uncoupled oscillator with self- 
feedback. In particular, considering only the longitu- 
dinal stabilization we have the limitation 5ft(Ao)T < 2. 
For details we refer to ref. [2^. However, there is still 
the transversal solution to be stabilized. The transver- 
sal equation coincides almost completely with the longi- 
tudinal one, but the delay term has the opposite sign. 
This means for our Floquet problem that, whenever the 
TT-torsion of a periodic orbit provides successful stabiliza- 
tion in the longitudinal direction, the inverted delay term 
and the torsion cancel each other in the transversal sub- 
space, just as if the orbit was torsion-free. The stability 
dilemma is exactly the same as for the fixed point of our 
system at half-period delay: For the presented coupling 
scheme at least one direction of the complete phase space 
remains unstable. But this is not an absolute restriction. 
Stable synchronized periodic oscillations may occur when 
the delay time is detuned from the orbit period, as can 
be seen in fig. 21 As mentioned in the preceding sec- 
tion, there could be an overlap between the domains of 
longitudinal and transversal stability, but only if the re- 
striction 5?(Ao)r < 0.326 is fulfilled together with a 7r/2 
torsion. For fixed points the overlapping domains were 
located around a quarter-period delay. For periodic or- 
bits the situation is more complicated, since the orbits 
also get deformed by the invasive feedback for t ^ 7i • T. 
But as a rough estimate the tt/2 argument holds for the 
stability domains of the period- 1 and period-2 orbits as 
well. 



VII. COMPLETE CHAOTIC 
SYNCHRONIZATION 

Complete synchronization of chaotic trajectories is ob- 
served for very small delay times r w but also for rather 



large delay times r « 3T, see fig. H] The first result 
is intuitively expected since it occurs close to the limit 
of conventional chaotic synchronization. This process is 
well established, and a slight change of the synchroniza- 
tion signal should not immediately spoil the whole mech- 
anism. The synchronization at t w 3T, however, was 
extremely surprising, even more when taking into consid- 
eration that the delay time exceeds the 'chaotic memory', 
i.e. the Lyapunov time, of the uncoupled chaotic oscilla- 
tor, by about factor 2. We evaluated this characteristic 
time to amount to 1.6T [2^. So it appears impossible 
that the tendency of nearby trajectories to exponentially 
diverge could be suppressed by coupling signals for such 
a large delay time [23]. What is also obvious for time- 
delayed chaotic synchronization, we do have invasive cou- 
plings. In fact, we do not synchronize the two original 
oscillators but - according to eq. ^ - two delay systems, 
which for increasing delay time and coupling strength will 
increasingly differ from the original oscillators. 

In order to analyze the observed synchronization phe- 
nomenon more in detail we again consider the tangent 
linear systems 

6u = T)i(C(t))Su + K(+Sur - Su) 

(11) 

which only differ from the previous forms, eqs. ^ and 
eqs. ([TU)) . by the common chaotic drive C{t)- The trajec- 
tory C{t) obeys the equation of motion 

C(t)=f(C(t))+K(C(t-r)-C(i)) , (12) 

which is the longitudinal part of eqs. ([5]) and includes the 
effect of the coupling parameters. Note that C{t — t) ^ 
C(i) for all T 7^ 0, which excludes any non- invasive de- 
layed couplings. The solutions (5u(t) and Sx(t) reflect the 
stability properties of C(i)- AH we can say is that there 
is some general tendency for bw{€) and 5v(t) to diverge 
(or converge) exponentially with the largest Lyapunov 
exponents A|| and A_l of the longitudinal and of the 
transversal subspace, respectively. A necessary criterion 
for chaotic synchronization is A_l < while A|| > 0. The 
challenging question is, can we present any limitations in 
terms of A|| and Ax for the domains of complete chaotic 
synchronization, as we could give for synchronized fixed 
points or periodic orbits? Certainly not in a strict way. 
Since the TLS, eqs. pT|) . are not driven by constant or 
periodic terms, which implies that Ju,- 7^ A • (5u and 
7^ A ■ 5v for any constant matrix A, we cannot di- 
rectly derive characteristic equations like eqs. ©. Yet, 
in the appendix, we develop an approximate form for a 
similar type of equations, which will allow us to specify 
limiting conditions for the domains of complete chaotic 
synchronization. 

A first guess of such a limiting condition can be ob- 
tained by means of a simple gedankenexperiment (see 
fig. [7]): Imagine a setup, which consists of two identi- 
cal oscillators like those used in our experiment. Let 
each of them include a self-feedback with fixed values of 



7 




FIG. 7. Gedankenexperiment: Two delay systems coupled 
via a non-invasive controller. 

time delay r and coupling K identical for both circuits. 
We consider this feedback to be part of each system, i.e. 
we have two identical delay systems, each of them char- 
acterized by their common largest Lyapunov exponent 
A|j. When A|| is positive, both delay systems are chaotic 
and, starting from neighboring initial conditions, they 
separate exponentially with this rate. The divergence 
of the trajectories can be compensated by an additional 
active or passive controller that applies an appropriate 
dissipative force to each oscillator. E.g. by means of a 
non- invasive control force K'(y — x) proportional to the 
difference between their instantaneous signals complete 
synchronization can be achieved, and the delay systems 
show the same dynamics as if they were not influenced 
by the controller. Now, let us assume that the controller 
has also a finite latency t' between measuring x and y 
and feeding the control forces back to the systems. Then 
the corresponding equations of motion take the form 

i = FK(x,x^)+K'(y^- -x^O 

y = FK(y,y.) + K'(x,, -y.O- ^ ' 

where Fk(x, x,-) and FK(y,yr) denote the delay sys- 
tems. It is a well-known fact in control engineering (see 
e.g. [2^ and references) that this latency r' has to re- 
main below some critical value in order to successfully 
stabilize the completely synchronized state x = y - i.e., 
returning to our previous notation, to stabilize the fixed 
point V = 0. If the instability were created by an ordi- 
nary system without delay one would expect that 

A t' < 1 (14) 

must be satisfied, where A is the real part of the corre- 
sponding exponent of the unstable fixed point the con- 
troller has to stabilize. By identifying the instability of 
the fixed point v = with the (longitudinal) Lyapunov 
exponent of our delay systems, A = A||, we obtain an 
estimate for the largest delay time r' accessible by non- 
invasive control. This analytical result is quite general 
and not restricted to the specific model of our experi- 
ment. 

If — like in our virtual setup — A itself is largely deter- 
mined by the intrinsic time delay r, which is of the same 
order as t', then we expect an interplay between the two 
different sources of delay, which might modify the con- 
ventional limitation (|14p . In particular, for r' = r and 
K' = K the eqs. (jl3|) just represent our original cou- 
pling scheme, eq. (jH). A more sophisticated approach 




FIG. 8. Numerical simulation of the product Ay ■ r shown 
in grayscale for synchronized solutions at t ~ ST. Hatched 
areas indicate absence of complete synchronization. 

presented in the appendix results in a moderate exten- 
sion of the simple estimate: 

A||r<1.28, (15) 

which means that the interplay between the non-invasive 
controller delay t' and the intrinsic delay r results in a 
30% increase of the stability limit. 

The idea of coupled delay systems provides the key for 
understanding the possibility of chaotic synchronization 
at large delay times, e.g. at r w 3T. Self-feedback in 
our delay systems creates a general tendency to reduce 
the grade of chaos and to make the systems more lami- 
nar. This is no surprise since the self-feedback acts just 
like the mechanism of Pyragas control. We have studied 
this reduction in detail by calculating numerical solutions 
C{t) of cq. and analyzing their longitudinal stabil- 
ity according to eq. pT|) . Results are presented in fig. [8] 
in grayscale. In order to limit our efforts we have only 
considered such sets of coupling parameters that result 
in complete synchronization. 

The synchronization patterns obtained by numerical 
simulations coincide almost exactly with our experimen- 
tal data, which means that the numerical model yields a 
quite realistic description of the experiment. Within the 
synchronization domains relation (|15p is clearly satisfied. 
For most of the parameter values, where an aperiodic 
synchronous signal was detected, the exponent is close 
to zero, so the trajectories are either weakly chaotic or 
quasiperiodic. The occurrence of quasiperiodicity might 
be related to the dynamics observed at the upper edge 
of a Pyragas control regime, where a controlled orbit be- 
comes unstable through a Hopf bifurcation and the dy- 
namic switches over to a torus. At those parameter sets, 
where the completely synchronized state remains clearly 
chaotic, the Lyapunov exponent A|| (fc, t) of this feedback- 
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disturbed chaos is still markedly smaller than the expo- 
nent of the original undisturbed oscillator, Aq. From our 
simulations Aq was evaluated to amount to O.lwo, which 
corresponds to a Lyapunov time Ta^ = I/Aq = 1.6 T. 
This means that for the fc — r area investigated the 
'chaotic memory' of the delay system has been extended 
in time by more than factor 5, which explains the possi- 
bility that complete chaotic synchronization can, in prin- 
ciple, occur for such large delay times r « 3T and for 
even larger ones. 

Finally we underline that complete chaotic synchro- 
nization, as observed for r w 3T, is not a singular event. 
For larger delay times we may observe even more ex- 
tended domains of complete synchronization - provided 
that the constraint ([15]) is still satisfied. In fig. ^ we 
present additional experimental data taken in the same 
setup, but for a different system parameter, a = 0.22. 
For this value of a the uncoupled oscillators arc still un- 
stable, but less chaotic than for our reference system 
(a = 0.265). The figure shows a lot of parameter do- 
mains in k and r where aperiodic synchronization was 
observed. In our automatic scans of the k—r plane we did 
not systematically distinguish between quasiperiodic and 
chaotic time series, but among those selected time series, 
which were analyzed in detail, we found many examples 
of complete chaotic synchronization. Complete aperiodic 
synchronization is again observed in a boomerang-shaped 
domain at r « 3T, but also in regions at higher multi- 
ples of T, where the shape of the domains differ qualita- 
tively for odd and even T. The increased size of the syn- 
chronization domains is obviously related to the weaker 
chaotic behavior of the uncoupled oscillators. Because of 
that the resulting longitudinal Lyapunov exponents A|| 
should be still smaller than for our reference system and 
- according to condition (jl5|) - synchronization should be 
accessible for even larger delay times. 




0.1 0.2 0.3 0.4 0.5 
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FIG. 9. (Color online) Dynamical states for system parameter 
a — 0.22. Labels as in fig. 2] 



VIII. CONCLUSIONS 

Two chaotic oscillators that are mutually coupled 
via time delayed signals show very complex patterns 
of completely synchronized dynamics, including fixed 
points, periodic orbits and quasi-periodic as well as 
chaotic trajectories. We have experimentally observed 
these synchronization phenomena in two delay-coupled 
Rossler-type electronic circuits and analyzed the syn- 
chronized states by means of analytical calculations and 
numerical simulations of the corresponding stability 
equations. We found that the equations determining the 
longitudinal and the transversal stability only differ by 
the sign of a single delay term, which, however, in the 
case of fixed points and periodic orbits results in limita- 
tions, which contradict or widely exclude each other. For 
instance, it is impossible to completely synchronize two 
Pyragas-controlled orbits, because the torsion, which 
favorites longitudinal stability, will prevent transversal 
stabilization and vice versa. For fixed points a similar 



phase argument excludes the simultaneous stabilization 
of longitudinal and transversal components for all but a 
very small set of coupling parameters. 

For chaotic trajectories such a dilemma does not exist, 
since merely transversal stability has to be achieved. 
Even more, if one of the stability exponents is allowed 
to be positive, this phase argument provides for a 
better chance to find the other one negative. Therefore, 
complete chaotic synchronization should occur more 
frequently, as long as the delay time does not exceed 
the limit of the chaotic instability. In fact, we observed 
complete synchronization of chaotic trajectories in our 
experiment for delay times far beyond the timescale of 
the uncoupled oscillators. This unexpected behavior 
could be attributed to the invasive character of our 
coupling scheme in the case of aperiodic trajectories. 
The coupling forces act like a virtual self-feedback, 
which has the tendency to reduce the grade of chaotic 
instability and to increase the Lyapunov time. As a 
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consequence complete chaotic synchronization could be 
achieved for very large delay times. We derived a general 
limitation for r connecting the transversal stability 
with the longitudinal Lyapunov exponent of the coupled 
chaotic systems. 

In spite of the simple and highly symmetric coupling 
scheme complete synchronization remains a rather ex- 
ceptional phenomenon, which shows a very complex de- 
pendence on the coupling parameters. The complexity 
seems to be characteristic for systems with time delay 
and might result from the interplay of diverging tenden- 
cies: to complicate the system by inflating the dimension 
of phase space, or to simplify the system through the ac- 
tion of an invasive self-feedback. We were able to present 
general arguments excluding or limiting complete syn- 
chronization for different types of dynamics. However, 
the detailed mechanism of transversal stabilization is not 
yet fully understood and has certainly to include more 
detailed information on the phase rotation (see [30l | ) and 
chaotic mixing of the components. Further investigations 
based on a stochastic modeling of the chaotic drive might 
shed more light on this problem. 



APPENDIX 

Starting from the tangent linear systems, eqs. ([TT|). 
we consider a small variation of the coupling strength 
k without changing the driving trajectory ^(Oj which is 
determined by cq. (|12p for some fixed value of k. Such 
a virtual step can be implemented only numerically, be- 
cause in experiment - due to the invasive character of the 
coupling - any variation of k would always influence C(t). 
By analogy with the procedure for periodic orbits in [27[ 
we introduce an approach for the Lyapunov exponents in 
the sense of a linear expansion in k around k = 



A = Ao + A: (x c 



-At 



(16) 



Here Aq depends only on the trajectory and is consid- 
ered as a hidden intrinsic property. If C were an unstable 
orbit of a single oscillator, Aq would correspond to the 
real part of its Floquet exponent in the absence of con- 
trol. Xi^y G 9're expansion coefficients quantifying the 
coupling efficiencies of time-delayed and instantaneous 
terms, respectively. In general wc have X which 
reflects the fact that for an aperiodic trajectory the de- 
layed coupling term cannot be mapped to the instanta- 
neous term, as mentioned above. Note that in contrast to 
the corresponding Floquet problem [23| eq. ([TB)) holds for 
Lyapunov exponents and - provided that k remains suffi- 
ciently small - all involved parameters are real. In terms 
of such an approach we estimate the limits of transversal 



stability for a given longitudinal instability A|| and delay 
time r. We start from eq. (fT6|) . which holds for both 
longitudinal and transversal dynamics 



A|| = Ao + fc(xe-^it"-r7) 



(17) 



Note that the symmetry of the longitudinal and transver- 
sal LTS results in identical expansion parameters x a-nd 77. 
The coefficient 7 includes a variation of sign and strength 
of the transversal coupling with respect to the longitu- 
dinal coupling and can also be applied to more general 
schemes than considered here [20[. The instantaneous 
terms in the coupling scheme can be easily included in 
the intrinsic Lyapunov exponent through the renormal- 



ization Aq — > Ag — 
transcendental eqs, 
W-function: 



r]k := Aq. Explicit solutions of these 
can be given in terms of the Lambert 



All = A' +-W(fcre-^n^ 

T V 

A_L = A^, + iw(7fcre-^«^ 



(18) 



The possibility of complete synchronization depends on 
the fact whether the global minimum of A_l can become 
negative. For given Aq the minimum is determined by the 
properties of the W-function and occurs at W{—c~^) = 
— 1, i.e., if k takes such a value that ^kTe~^°'^ — —1/e 
we have 



A, 



and correspondingly. 



An - A' 



T 



1 / -1 

- W — 

T V ^7 



(19) 



(20) 



Eliminating Ag from both eqs. we finally arrive at the 
relation 



A±,min — A|| 



1 



1 + W 



-1 

e7 



(21) 



which connects the minimum value of A^ with the corre- 
sponding longitudinal exponent A|| and the delay time t. 
Stable transversal dynamics only occurs for A^^min < 0, 
which directly leads to a general limitation for the delay 
time: 



All • r < 1 -f W ( — 

* 67 



(22) 



In our present coupling scheme we have 7 = — 1, which 
results in the limiting condition A|| ■ t < l + W(+l/e) « 
1.28. 
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